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Abstract 

A qqq BSE formalism based on an input 4-fermion Lagrangian of 'current' 
u, d quarks interacting pairwise via a gluon-exchange-like propagator in its non- 
perturbative regime, is employed for the construction of a relativistic qqq- wave func- 
tion under the Covariant Instantaneity Ansatz (CIA). The chiral invariance of the 
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input Lagrangian is automatically ensured by the vector character of the gluonic 
propagator, while the 'constituent' masses are the low momentum limits of the dy- 
namical mass function m(p) generated by the standard mechanism of DB\S in the 
solution of the Schwinger Dyson Equation (SDE). The CIA gives an exact reduction 
of the BSE to a 3D form which is appropriate for baryon spectroscopy, while the 
reconstructed 4D form identifies the hadron quark vertex function as the key ingre- 
dient for evaluating transition amplitudes via quark- loop integrals. In this paper 
the second stage of this 'two-tier' BSE formalism is extended from the 4D gg-meson 
to the 4D qqq-haiyon vertex reconstruction through a reversal of steps offered by 
the CIA structure. As a first application of this 4D qqq wave function, we evaluate 
the quark loop integrals for the neutron (n) - proton (p) mass difference which re- 
ceives contributions from two sources : i) the strong SU(2) effect arising from the 
u — d mass difference (4 MeV); ii) the e.m. effect of the respective quark charges. 
The resultant n — p difference works out at 1.28 MeV (vs. 1.29 expt), with two free 
parameters Cq,u}q characterizing the infrared structure of the gluonic, which have 
been precalibrated from a common fit to qq and qqq spectra as well as several other 
observable quark loop integrals. A formal derivation, based on Green's function 
techniques for 3 spinless quarks, of the CIA structure of the 4D qqq-b&ryon vertex 
function as employed in the text, is given for completeness in Appendix B. 
PACS : 11.10 st ; 12.35 Ht ; 12.70 + q 
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1 Introduction: Unified BS Dynamics of 2- and 3- 
Quark Hadrons With 3D Kernel Support 

Soon after the advent of the Faddeev theory [1], the relativistic 3-body problem [2] at- 
tracted instant attention as a non-trivial dynamical problem, as distinct from earlier 
"kinematical" attempts [3] at a relativistic formulation of its wave function. In this re- 
spect the relativistic 3-baryon problem had been more of academic than practical interest 
(until the 'pion' got involved as a key ingredient), but the situation changed qualitatively 
when this 3-body problem started being viewed at the quark level. Looking back after 
25 years it appears that the first serious attempt in this direction was made by Feynman 
et al [4] who gave a unified formulation of both the qq (meson) and qqq (baryon) prob- 
lems under a common dynamical framework, bringing out rather sharply an underlying 
duality between these two systems which in turn signifies a more basic duality between 
a qq diquark [5] and a q antiquark. Indeed the diquark description is quite compact and 
adequate for many practical purposes involving the baryon, but the more microscopic 
qqq description which brings out the fuller permutation (S3) symmetry in the baryon is 
necessary for the actual details of a full-fledged dynamical treatment [4]. 

Although the FKR theory [4] marked the first step in this direction, it suffered from 
an inadequate treatment of the time-like d.o.f. which showed up in several ways. The 
latter has by itself a long history in terms of attempts at formulating the Bethe-Salpeter 
Equation (BSE) for qq qqq systems throughout the Seventies under the Instantaneous 
Approximation (IA), as has been reviewed elsewhere [6]. The normalization of the BS 
wave function (gaussian) has been given by Tomozawa [7] under some special assumptions 
which however are not general enough to be adapted to any broader dynamical BSE 
framework going beyond the IA. 

For several years we have been involved with a certain formulation of BS-dynamics 
for both qq and qqq systems within a common unified framework (to emphasize their 
underlying duality), designed to address their spectroscopy on the one hand, and a self- 
consistent treatment of various quark-loop amplitudes in terms of their respective hadron- 
quark vertex functions on the other. The "spectroscopy" aspects are addressed through 
the 3D reductions of the qq BSE [8] and the qqq BSE [9], to compare with observed 0(3) 
spectra [10], while the loop-integral aspects of transition amplitudes show up through the 
reconstructed vertex functions of the 4D BSE [11]. An exact interconnection between the 
two forms was achieved through the ansatz of a 3D support defined covariantly in the BS 
kernels for the qq [12] and qqq [13] systems. 

The "Covariant Instaneity Ansatz" (CIA) [12] which has been at the root of this 
"two-tier" philosophy, and is also supported by other considerations [14] based on the 
Markov- Yukawa transversality condition [15], seems to offer a possible Lorentz covariant 
way to reconcile the apparently conflicting demands of 3D spectroscopy [10] with the 4D 
structure of quark-loop amplitudes. The effectiveness of the CIA in giving a concrete 
shape to such a "two-tier" philosophy of spectra-cum-loop integrals was summarised in a 
semi-review [13] in the form of appropriate BSE's for qq and qqq systems with vector-type 
kernels [6] with 3D support, albeit with slight modifications [6] in the respective BSE 
structures to facilitate greater 'manoeuvreability', in the spirit of similar efforts [16] in 
the past. Further, the observed spectroscopy [10] is well satisfied on both qq [8] and qqq 
[9] sectors with a common set of parameters for the respective kernels (the qq kernel has 
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just half the strength of the qq kernel due to color effects), so that the respective vertex 
functions are entirely determined within this formalism. 

The other aspect of this 'two-tier' formalism concerns the crucial property of chiral 
symmetry and its dynamical breaking. The first part (chiral symmetry) is ensured without 
extra charge by the vector character of the kernel that had been present all along in this 
program [6,11], since the BS-kernel is a direct reflection of an effective 4-fermion term in 
the input Lagrangian. Indeed the vector type character of the latter lends a natural gluon 
exchange flavour to such a pairwise interaction among 'current' (almost massless) u, d 
quarks at the Lagrangian level. This structure is quite general [17], and can be adapted 
to the QCD requirements on the gluonic propagator involved in the pairwise interaction 
kernel. Of this, the perturbative part (which is well understood) is quite explicit, but the 
non-perturbative (infrared) part is not yet derivable from formal QCD premises. It can 
nevertheless be simulated in a sufficiently realistic manner at the phenomenological level 
[6,19], so as to satisfy the standard constraints of confinement as well as explicit QCD 
features [20] in terms of a basically 3D BSE kernel structure. 

The second part, viz., dynamical breaking of chiral symmetry {DB\S) is implemented 
via the Nambu-Jonalasino mechanism [21] whose full-fledged form amounts to adopting 
the 'non-trivial' solution of the Schwinger- Dyson Equation (SDE) derived from a given 
input, chirally symmetric Lagrangian with current quarks. A mass function m(p) [17,18] 
is thus generated whose low-momentum value may be identified with the bulk of the 
'constituent' mass (m q ) of the u, d quarks. This accords with Politzer additivity [22], viz., 
m q = m(0) + m c ; where m c , the current mass, is small. This was also shown in the context 
of a BSE-cum-SDE treatment [23] within the CIA formalism [12]. Thus formally the BS- 
kernel may be regarded as a non-perturbative gluon propagator [23] in a BSE framework 
involving the dynamical/constituent mass [19, 21-24] in the quark propagator. 

To recapitulate, the CIA which gives an exact interconnection between the 3D and 4D 
forms of the BSE, provides a unified view of 2- and 3-quark hadrons, its 3D reduction being 
meant for spectroscopy [8-9], and the reconstructed 4D form [12,13] for identifying the 
respective hadron-quark vertex functions as the key ingredients for evaluating 4D quark- 
loop integrals. The formalism stems from a strongly QCD-motivated Lagrangian with 
current quarks whose pairwise interaction is mediated by a gluonic propagator in its non- 
perturbative regime. The QCD feature of chiral symmetry is ensured by the vector nature 
of this interaction, while its dynamical breaking is the result of a non-trivial solution of the 
SDE [17,23]. Thus, unlike in conventional potential models [25], the constituent mass so 
generated is not a phenomenological artefact, but the result of a self-consistent solution of 
the SDE [17, 19, 23], wherein the (constituent) mass normally employed for spectroscopy 
[8-9] matches with the output dynamical mass at low momentum [23]. Thus there are only 
two genuine input parameters Co, u>o, that characterize the (phenomenological) structure 
of the non-perturbative gluon propagator which serves for both the 2- and 3-quark spectra 
in a unified fashion [8-9]. In this formalism, these two constants play a role somewhat 
similar to that of the (input) 'condensates' in the theory of QCD sum rules [26]. 

Before proceeding further, let us pause to compare this approach with other dynamical 
methods, e.g., chiral perturbation theory [27] which has more explicit QCD features, albeit 
in the perturbative regime, leading to expansions in the momenta. This is a powerful 
theoretical approach employing the (chiral) symmetry of QCD; its essential parameters 
are the current quark masses, and the method works very efficiently where its premises 
are logically applicable. Thus it predicts the ground state spectra of light quark hadrons, 
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including their mass splittings due to strong and e.m. breaking of SU(2), but not the 
spectra of L-excited hadrons. The latter on the other hand demand a "closed form" 
approach to incorporate the "soft" off-shell effects which in turn require a non-trivial 
handle on the infrared (non-perturbative) part of the gluonic propagator, something which 
the present state of the QCD art does not yet provide. Thus one needs a phenomenological 
input even in standard BSE-SDE approaches [18], as discussed elsewhere [23]. The chiral 
perturbation theory [27] also lacks this vital ingredient, as seen from the absence of form 
factors in its 'point' Lagrangians [27] with at most derivative terms. This shows up, e.g., 
through its inability to predict L-excited spectra, and finer aspects (such as convergence) 
of 4D quark-loop integrals which depend crucially on these "off-shell" features. Physically 
this amounts to the absence of a 'confinement scale' which governs these form factors. In 
other BSE-cum-SDE approaches [17-19], including the present 'two-tier' CIA formalism 
[12,23], this 'scale' is an integral part of the structure of the non-perturbative part of 
the gluon propagator [19,23], with a built-in QCD feature of chiral symmetry and its 
dynamical breaking through the non-trivial solution of the SDE [17,19,23]. This not only 
facilitates the prediction of L-excited spectra [19, 8-9] but also provides a form factor for 
the hadron-quark vertex function which greatly enhances its applicability to various 4D 
quark- loop integrals; see [12, 23-24, 28]. 

After this clarification on the philosophy of this two-tier BSE approach, vis-a-vis some 
others [26,27], we may now state the objective of the present paper: A typical application 
of the 4D baryon-ggg wave function reconstructed [13] from the 3D qqq BSE, as a 3- 
body generalization of the corresponding gg-meson problem [12]. Unlike the 2-body case, 
however, where the steps are exactly reversible [12], such reconstruction in the 3-body 
case involves a loss of information on the 4D Hilbert space for a 3-body system, so that 
the reversal of steps is in principle not unique, and requires a ID 5-function to fill up 
the information gap between 3D and 4D Hilbert space which may be directly attributed 
to the CIA ansatz of a 3D support to the pairwise kernel. The 2-body case just escapes 
this pathology as it represents a sort of degenerate situation, but the price of a 3D kernel 
support must show up in a reconstruction of the 4D BSE from its reduced 3D form in 
any (n > 2)-body problem [29]. A plausible 'CIA' structure for the 4D qqq wave function 
was suggested in [13] in a semi-intuitive fashion, but a more formal mathematical basis 
has since been found [29] through the use of Green's function techniques, so that the 
reconstructed 4D form reduces exactly to the (known) 3D form as a consistency check 
[29]. The final result, which is almost the same as the earlier conjecture, eq.(5.15) of [13], 
except for a constant that does not affect the normalization, contains a ID (^-function 
corresponding to the on-shell propagation of the spectator between two successive vertex 
points. As explained in detail in [29], this ID 5-function must not be confused with 
any signature of "non-connectedness" in the 3-body wave function [30], since the 3D 
form is fully connected. A better analogy is to the 'scattering length approximation' 
to the n — p interaction, characterized by the appearance of a (Fermi-type) 5-function 
potential, in estimating the effect of chemical binding on the scattering of neutrons by 
a hydrogen molecule [31]. In any case the ID 5-function appearing in this structure is 
entirely innocuous since it gets integrated out in any physical (quark loop) amplitude 
including the BS normalization (see Sec. 2 below). 

The application chosen for the baryon-ggg vertex function is to the n-p mass difference, 
on closely parallel lines to the meson case [28]. To racall the physics of the n-p mass 
difference, this quantity receives contributions of opposite signs from two distinct sources 
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: i) a positive one from the strong SU(2) d — u mass difference; 
ii) a negative one from e.m. splittings. 

For the sake of completeness, Appendix B of this paper gives the main steps of the 
derivation [29] for the 4D structure, eq.(5.15) of [13], of the baryon-ggg vertex function by 
the Green's function method for 3 spinless quarks, to be employed in this paper. In Sec. 2 
we collect the various pieces of this quantity with the inclusion of the spin and isospin 
d.o.f., on the lines of [32]. Thus equipped, we outline the main steps leading to an explicit 
evaluation of the normalization integral, using Feynman diagrams shown in figs.l(a.b,c). 
A complex basis [9, 33, 34] for 3D momentum variables facilitates the evaluation of the 
resulting 3D x 3D integrals, after the time-like momenta have been eliminated by 'pole' 
integrations on identical lines to the corresponding qq problem [12,23,24]. In Sec. 3 we 
evaluate the 'shift' in the nucleon mass due to strong SU(2) breaking, by inserting a 
mass shift operator —SmT 3 ^/2 in place of i^^Ci at each of the corresponding 7- vertices 
of figs.l( shown in figs. 2.(a,b,c). Here 8m = 4 MeV is the 'standard' d-u 

mass difference [24,28] taken as the basic input. The condensate contribution is neglected 
as it was found to be negligible from a similar calculation of the SU(2) mass splittings 
in pseudoscalar mesons [28]. Sec. 4 sketches the evaluation of the e.m. contribution in 
accordance with the diagrams of fig.3(a,b,c), while the details of the approximations 
employed are collected in Appendix A. Sec. 5 summarises our findings and conclusions 
vis-a-vis other methods. 

2 Normalization of the Baryon-qqq Vertex Function 

To outline the structure of the baryon-ggg vertex function from a CIA-governed BSE 
[12-13], we shall generally follow the notation, normalization and phase convention for 
the various symbols as given in [13], but adapted to the equal mass kinematics (m x = 
™>2 = 1^3 = Tn q ). The SU(2) mass difference 5m AMeV) between d and u quarks will 
be taken into account only through a 2-point vertex \—5mT^/2] inserted in the quark 
propagators in figs. 2 (in place of «7 M ej for a photon), but not in the structure of the vertex 
function. The vertex function is written in three pieces in each of which one quark plays 
the role of the 'spectator' by turn. For the spin structure (not given in [13]) we employ 
the convention of [3] which was extended in [32] to incorporate the ^-symmetry for the 
spin-cum- isospin structure in the Verde [35] notation [36] . The full 4D BS wave function 
^ reads as [13,32,33] : 

^AxA 2 A 3 = (v 1 + r 2 + r 3 ) x [ X V + X n r}/V2; (2.1) 

Ai = m q 2 + Pl 2 ; (i = 1,2,3). (2.2) 

Here x' an d x" are the relativistic "spin" wave functions in a 2-component mixed sym- 
metric S3 basis which for a 56 baryon go with the associated isospin functions 0' and (/>" 
respectively. These are given by [3,17] : 

[x'Wa = [(M - i r P)i l5 C/y/2]^ x U(P)J(2M) (2.3) 

[ X "K Q = [(M - i r P) llt C/y/6]^ x iw„U(P)J(2M) (2.4) 

in a spinorial basis [3,32] in which the index a refers to the 'active' quark (interacting 
with an external photon line, fig.l), while /3, 7 characterize the other two, with the further 
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convention that 7 refers to the "spectator" in a given diagram, fig.(l). The 'hat' on 
7 signifies its perpendicularity to P M , viz., 7.P = 0. The notations in eqs. (2.3-4) are 
standard, with a common Dirac basis for the entire structure, and 'C is the charge 
conjugation operator for quark #3 in a 23-grouping [3,32]. P M is the baryon 4-momentum, 
U(P) is its spinor representation, and (M — i r y.P)/(2M) its energy projection operator 
[3,32]. Further, because of the full S^-symmetry of the last factor in (2.1), the (1,2,3) 
indices can be permuted as needed for the diagram on hand. Thus in fig. la, #l(a) 
interacts with the photon ; #2(/3) is the quark which has had a 'last' gg-interaction with 
#l(a) before emerging from the hadronic 'blob', while #3(7) is the spectator [32]. In 
fig. lb, the roles of #1 and j^l are reversed so that, of the two 'active quarks' #1 and 
#2, #2(a) now interacts with the photon, #l(/3) has had the last qq- interaction with 
#2(a), while #3(7) still remains the 'spectator'. These roles are cyclically permuted, 
with two more such pairs of diagrams, fig.lc), to give an identical chance to each of the 
quarks in turn [32]. Thus there are 3 such pairs of diagrams, of which only one pair is 
shown. An identical consideration applies to figs.2(a,b) with i^^ei replaced by {—8mT 3 /2) 
consistently. The spatial vertex functions Tj are given for i = 3 by [13] : 

T 3 = N B [D 12 <t>/2m] x ^/[27r5(A 3 ).A 3 ] (2.5) 

where is the full, connected qqq wave function in 3D form, and D 12 is the 3D denominator 
function of the (12) subsystem . The second factor represents the effect of the spectator 
[13] whose inverse propagator Pp _1 (p 3 ) off the mass shell is just A 3 , eq.(2.2). The 
main steps leading to this unorthodox structure which has been derived recently via 
the techniques of Green's functions [29], are sketched for completeness in Appendix B. 
As already noted in Sec.l, and again explained in Appendix B, its peculiar singularity 
structure in the form of a "square-root" of a ID 5-function stems from the CIA ansatz 
of a 3D support to the pairwise interaction kernel, but it is quite harmless as the former 
will appear in a linear form in the transition amplitude corresponding to any Feynman 
diagram as in figs. 1-2. The complete expressions for D 12 and are given for the equal 



mass case (with #3 as spectator) by [13] (see also [9]): 

D 12 = A 12 (M-cu 3 ); A12 = 2uj 12 2 - M 2 (l - z/ 3 ) 2 /2 (2.6) 

^i2 2 = m q 2 + qf 2 ; 2qf 2 = p£ - p$ (2.7) 

= e -(pi+Pl+pl)l^ 2 = e -p/W 2 (2.8) 

(see further below for the definition of p). 

p»= P S+ Pl .PP,/M 2 ; #+$+$ = (2.9) 

uj 2 = m 2 + p 2 ; u 3 = uj 3 /M(onshell) (2.10) 
The /5-parameter is defined sequentially by [8,9]: 

p 4 = ^Mcu 2 a s (l -m q /M) 2 {M- < to >); < u > 2 = m 2 + 3/3 2 /8 (2.11) 
9 

a s =a s -2MC M _ < - > , (2.12) 
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^ = 291n (M - < " >) ; (2.13) 

«s J^QCD 

A QCD = 200MeV; cu = l58MeV; C = 0.29 (2.14) 

The normalization Nb, eq.(2.5), is given in accordance with the Feynman diagrams l(a,b) 
by the 4D integral ( c.f.[32]) : 

123 " 2A 2 A 3 

+ (1^2) (2.15) 



iPjM = Y,j d ^ d ^^T [< ^WW* > +3 < 0"I(23)\a(1)'WI0" >] 



where the matrix element for fig. la is organized as a product of two spin-factors : a 
'23-element' expressed as a Dirac trace over the indices f3, 7 ; and a '1-element' (with 
suppressed index a). The associated isospin functions are shown according to (2.1). 
The contribution of fig. lb is shown symbolically by 1 <^ 2, while ^123 indicates the sum 
over all the 3 pairs cyclically. In representing eq.(2.12) we have dropped 'cross-terms' 
like T*Tj, where i 7^ j, since the presence of a ^-function in each Tj ensures that a 
simultaneous 'on-shell' energy conservation of i 7^ j spectators is not possible [32]. The 
various pieces of the matrix elements in (2.14) which can be read off from fig. la in terms 
of the spin functions (2.3-4) are as follows: 

(1)'W = U(P)S F (p l )t lfl e l S F (p 1 )U(P) (2.16) 

iS F ^ 1 (p) = m q + i>y.p (2.17) 

(1)'W = U(P)t% l5 S F (p 1 )t^e 1 S F (p 1 )t 1 ^ x U(P); (2.18) 

(23)' = Tr[C- 1 75 (M - i 7 .P)(m, - i 7 .p 2 )(M - ?7 .P) 75 (m g + ? 7 .p 3 )C]/8M 2 (2.19) 

(23)\ A = Tr[C^%(M - i r P)(m q - i 7 .p 2 )(M - i~f.P)%(m q + t 1 .p 3 )C]/8M 2 (2.20) 

The 'strength' of the (zero- momentum) 'photon' coupling to the quark line pi can be 
chosen in several ways [37]. We take here the simplest possibility, viz., ei = 1/3 each. 
The isospin matrix element is first eliminated according to [38]: 

< 0'|l|0' >=< 0"|1|0" >= 1 (2.21) 

< 0'|r 3 (1 V >= -3 < 0"r 3 «|0" >=< r 3 > M (2.22) 

Eq.(2.20) suffices for (2.14), while (2.21) will be needed for the u-d mass difference operator 
—SmTs^/2 ; see Sec. 3. Next, the evaluation of the traces in (2.15-18) is straightforward, 
after noting that (2.15-16), after spin-averaging, are expressible as traces. The results are 

(23)'6 uX = (23)\ A = (m q + Mu 2 )(m q + Mu 3 )6 vX (2.23) 

(1)^, A = (1)\ A;M = \2Mv x {m q + Mv x ) + A 1 ]^? /1 /(MA 1 2 ) (2.24) 
where 6 is a covariant Kronecker delta w.r.t. P M , viz., 

e vX = 9 vX = 5 uX - P U P X /P 2 ; (P 2 = —M 2 ) (2.25) 

Collecting all these results and simplifying we get 

Nb~ 2 = £ / d% {m % +UJ3) x / d% 2 D 12 2 <t> 2 [ ei h + e 2 / 2 ] (2.26) 
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2mh = J Mda X2 \2Mv x (m q + Mv x ) + Ai]/(M Ax 2 A 2 ) (2.27) 

where we have "cashed" the <5(A 3 )-function arising from |r 3 | 2 against the time-like com- 
ponent of d A p 3 , and used the results 

A12 = d 3 q 12 Mda 12 ; i/ 1>2 = (1 - u 3 ) ± a 12 (2.28) 

The integration over da 12 involves single and double poles arising from the propagators 
Ai j2 _1 in (2.26), while the value of v 3 is taken 'on-shell' at co 3 /M after the 5 (A3)- function 
has been cashed. The result of a basic ^-integration is 

J Mda 12 A l - 1 A 2 ~ 1 = 2m/D l2 (2.29) 

from which others can be deduced by differentiation under unequal mass kinematics, or 
directly through a 'double pole' integration. The net result for l x + 1 2 , eq.(2.26), is given 
in eq.(2.41) below. Further, the individual terms of the summation ^123 i n (2.25) are 
fixed by the values chosen for (which need not be specified in advance, as they can be 
adapted to other conventions too [37]; see Sec. 3). 

The integration in (2.25) can be considerably simplified in a complex basis [9,33] 
defined (in momentum space) by : 

V2 Zi = & + im; V2z* = & - Vi ; (2.30) 

= pu - P2U 3^ = -2p 3i + pu + p 2i ; (2.31) 

where we now employ the alternative notation pu for pi, in view of its basically 3D 
content. In terms of Zi and z*, the 6D integration in (2.25) is expressed as 

d 3 p 3 d 3 q 12 = (V3/2) 3 d 3 £d 3 r] = d 3 zd 3 z* (2.32) 

The further representation [9,33] 

d 3 zd 3 z* = (dz + dz_*)(dz-.dz + *)(dz 3 dz3*) (2.33) 

where 



V2z + 


= Ri e Wl - 


V2z-* 


= R l e~ i91 


(2.34) 


V2z_ 


= R 2 e if?2 ; 


V2z + * 


= R 2 e~ id2 


(2.35) 


V2z 3 


= R 3 e 103 ] 


V2z 3 * 


= R 3 e~ i93 


(2.36) 



reduces the 6D integration (2.32) merely to ix 3 dRi 2 dR 2 2 dR 3 2 , since the 6i -variables (not 
Euler angles!) are not involved in the integrands encountered, and just sum up to (27r) 3 . 
The positive variables Ri, (i = 1,2,3), are related to the £,i,r]i variables by 

p = Rl + R\ + R\ = £ 2 + rf = 2 Zi z* (2.37) 

To convert the variables Ui that appear in the integrals (2.28) in terms of the #1,2,3 
variables is a straightforward but tedious process which can be somewhat simplified in 
terms of the intermediate variables £ 2 — rj 2 and 2C,.i] which form a [2,1] representation [35] 
of ^-symmetry at the 'quadratic' level. Now because of the full S3 -symmetry of the 6D 
integral (2.32), together with the (fortunate) circumstance of equal mass quarks in the 
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problem on hand, the integrand as a whole is S3 -symmetric which permits the following 
simplification: Each of the quantities p 2 and q 2 inside (2.32) can be expanded as 

Fi,2 = P/2 + (e - V 2 )/* ± Pi = p/2 - (e - V 2 )/2 (2.38) 

f 12 = 3e/4 = p/2 + (e-V 2 )/4 (2.39) 

In all these terms the principal quantity is p/2, while the resultant effects of the mixed- 
symmetric corrections will show up only in the fourth order, etc. In the present case of 
equal mass kinematics it is a good approximation to neglect the latter terms, as has also 
been found for the qqq mass spectral results [9], so that all quantities are expressed in 
terms of p only : 

^1,2,3 ~ W12 ~ u P \ u p 2 = m 2 q + p/2 (2.40) 

D 12 » 2(M - co p )[uj 2 p — (M — uj p ) 2 /4]. (2.41) 

The rest of the integration is expressed entirely in terms of the p-variable, with the 
resultant 6D measure given by 

J d 3 p 3 d 3 q 12 F(p) = (ttv / 3/2) 3 j p 2 dp/2F(p) (2.42) 

These considerations suffice for evaluating the integrals Ii and 7 2 whose resultant value 
is now given for = 1/3 by : 

D 2 12 (h+I 2 ) = [m 2 q +m q (M-uj p ) + (M-uj p ) 2 /4}x(M-co p ) 3 /u p +D 12 (2m q +M-u; p ) (2.43) 

Substitution in (2.25) yields N B directly. The numerical values are given collectvely at 
the end of Sec. 4. 



3 Strong SU(2) Mass Difference for the Nucleon 

This calculation is on almost identical lines to Sec. 2, except for the substitution ie 1 / ~f fl 
to —5mT3^/2 in figs.l(a,b) to give figs.2(a,b) which represent the effect of insertion of 
a 2-point vertex in a quark line. Indeed we can directly start from the counterpart of 
eq.(2.15) which gives the 'strong' mass shift as: 

** M « = E/ d ^ 2d ^oWk X [< 0,|(23) ' (1) ' 10, >+ l < ^"K 23 )"^(1)",aI0" >] 

+ (1^2) (3.1) 

where we have now employed eq.(2.21) for the isospin factors, and the counterparts of 
(2.16) and (2.17) are respectively 

(1)' = f/(P)^(p 1 )[-W 3 (1) /2] 1 S F (p 1 )f/(P); (3.2) 

(1)\ A = f/(P)^75^(pi)[-W3 (1) /2]5 F ( Pl )2757Af/(P) (3.3) 

while the definitions (2.18) and (2.19) remain unaltered. As a result, eq.(2.22) remains 
valid, while the counterpart of (2.23) becomes 

(1)'0„ A = -3(1)\ A = [2m q (m q + Mv x ) - A 1 ]e uX (-5m/2)/A 2 1 (3.4) 
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Carrying out the o?<7i 2 -integration, the result for 5M st is now given by the counterpart of 
(2.25), viz., 

5M st = 3N 2 J d% x J d^D^^J, + J 2 ](-W 3 /6) (3.5) 

in the form of an isospin operator "r 3 " for the nucleon, where we have represented the 
effect of by a factor of "3" , and 

£> 2 2 [Ji + J 2 ] = ^[{m q + -{M-u p )){M-uu p )[2uu 2 + m q {M-uj p )} 

+ (M - uOpY\m q + (M — uOp)/2} 2 + A 12 (m 2 - cu 2 p ) 
+(m - LO p ){m q + (M — cj p )/2)(cu 2 + m g (M - co p )/2) 
+A 2 2 (l-m 9 (M-^)/^)/2] (3.6) 

as the exact counterpart of (2.42) under the same approximation. It is seen from (3.5) 
that the difference (n — p) is positive. 

4 E.M. Mass Difference for the Nucleon 

The diagrams for the e.m. mass difference are given by figs. 3 (I, II, III) for a proton (uud) 
configuration to illustrate the underlying topology in accordance with the roles of the 
'active' and 'spectator' quarks in turn, as explained in Sec. 2. In each of these diagrams, 
two internal quark lines are joined by a photon line. The e.m. vertex at quark j^i has the 
strength e[l + 3r 3 ^]/6 from which the isospin matrix elements of a product of two such 
factors (shown for fig.3.III) have the forms 

<0';0"|(l + 3r 3 (1) )/6x (l + 3r 3 (1) )/6|0';0" > (4.1) 

for the proton (uud) cnfiguration shown in III with #3 as spectator, but in a basis (1;23) 
(which is consistent with the spin basis, eqs.(2.3-4)), corresponding to fig. la, viz. [36,38]: 

\<f)' >= ui(u 2 d 3 - u 3 d 2 ) / V2; |0" >= (-2d 1 u 2 u 2 , + u 1 d 2 u 2 , + u^d^/VS (4.2) 

We note in parentheses that in fig. 3. Ill, the interchange of the two 'active' quarks #1 and 
#2 does not give a new configuration, unlike in figs.l and 2; ((a) versus (b) configurations). 

It is now easy to check that the matrix elements <>' and <> " of (4.1) are 1/9 and 
—1/9 for the proton configuration. After doing the corresponding neutron case, the two 
results may be combined in the single operator forms [38]: 

< . >'= (1 + 3r 3 )/36; < . > " = (1 - 5r 3 )/36 (4.3) 

where r 3 is the isospin operator for the nucleon as a whole [see eq.(3.5)], to be sandwiched 
between the neutron and proton states. The resultant isospin factor is then 

e 2 [< .>' + <.> "]/2 = e 2 (l - r 3 )/36 -e 2 r 3 /36 (4.4) 

After this book-keeping on the charge factors we can drop the isospin d.o.f. \<p > from 
the qqq wave function and, on the basis of the equality of the (.)' and (.)" contributions 
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(2.22-23) for the spin matrix elements, it is enough to work with the (.)' type to represent 
the full effect. Collecting these details, the net isospin contribution to the e.m. (n — p) 
mass difference is just e 2 /18, which (of course) comes out with the correct (negative) sign 
in the resultant e.m.contribution to the total n — p difference after all the phase factors in 
the orbital-cum-spin space have been taken into account. The complete e.m. self energy 
of the nucleon (with operator r 3 ), with fig.3.III as the prototype, is now given by 

6NT = ]r[-e 2 T 3 /36]/(27r) 4 / d% d% 2 d% 2 T;T' 3 k- 2 x [23];[1];/(A 3 A 2 A' 2 ) (4.5) 

123 

where the various momentum symbols are as shown in fig. 3, with the primed quantities 
referring to the vertex on the right, but otherwise written in the same convention as in 
eqs.(2.5-10). The symbols within square brackets are analogous to (2.16-19): 

[lr^uiP^sM^s^uipy, (p' = p) (4.6) 

Tt 

[23]; = ^pi C ^(M - i 7 .P')K - z7-P 2 )*7mK " *7-P 2 )(M - i r P)-y 5 (m q + H-Ps)C} 

(4.7) 

And the product of (4.6) and (4.7) works out as 

ME = (m ^ 3) [-(Ai + A; - k 2 )(A 2 + A' 2 - fc 2 )/4 

-(Ai + A; - k 2 ){m q u 2 + m q J 2 + 2uj 2 uj' 2 )/2 
-(A 2 + A' 2 - k 2 )(m q uji + m q J x + 2u 1 J l )/2 

-(Ai + A 2 ){m q + u[)(m q + cu' 2 )/2 - (A[ + A' 2 )(m q + ui){m q + u 2 )/2 
-(Ai + A 2 )(m 9 + 0Ji){m q + u' 2 )/2 - (Ai + A' 2 )(m q + ^)(m g + u 2 )/2 

+ ( m ? + _ P3 ~ /C ) 2 )[( m ? + ^l)( m 9 + ^ 2 ) + ( m Q + ^l)( m 9 + ^2) 

+ (m q + Ui)(m q + lu 2 ) + (m q + ^)( m g + ^2)]] ( 4 - 8 ) 

Some features of this "master" expression may be noted. There is a 'natural factor- 
ization'in the variables q\ 2 and q' l2 , except for the photon propagator k~ 2 , (k = q\ 2 — q[ 2 )- 
Further, the two blobs are connected together by the 'spectator 'variable p 3 which is on 
the mass shell due to the presence of TgTg in eq.(4.3). 

The time-like (pole) integrations over each of d<j\ 2 and a' 12 can be carried out exactly a 
la (2.28) and its derivatives, since the 3D vertex function Di 2 <p in T 3 does not involve o"i 2 , 
etc. After this step q± 2 , q' 12 and p 3 are the 'right' 3D variables for the 'triple integration' 
whose essential logic may be stated as follows. The main strategy is to decouple the q 
and q' variables from the photon propagator k through the following device [19] : 

Since k is basically space- like, it is a good approximation to replace k~ 2 by k 2 which 
equals (qi 2 — q' 12 ) 2 , and drop the angular correlation in the two q- momenta (since the 
error in this neglect is zero in the first order [28]). Next we use the inequality [28] 

(a 2 + b 2 )- 1 < (2ab)- 1 ; a -> \q 12 \,etc (4.9) 

which ensures the necessary factorizability in the q-variables.In principle the corrections 
to this inequality can be calculated since the neglected term is approximately equal to 
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— (a — 6) 2 /(4a 2 6 2 ) which is still factorizable, but this refinement is unnecessary in view of 
the smallness of the e.m. effect itself. After this simplification the rest of the integration 
procedure is straightforward since the q and q' integrations can be done analytically, 
and only a ID integration over |p 3 | remains for numerical evaluation. The necessary 
expressions are collected in Appendix A and the numerical results for all contributions 
are given as under. 

The key parameters are the quark mass m q and the size parameter j3 2 , the latter being 
determined dynamically through the chain of eqs.(2.11-14). As noted in Sec.l already, 
the mass m q which is usually called the 'constituent' mass, should be viewed as the sum 
of the (flavour independent) 'mass function' m(p) for small p, plus a small "current mass" 
m c , in the spirit of Politzer additivity [22]. The mass function m(p) was generated in this 
BSE-cum-SDE framework through a Dynamical Chiral Symmetry Breaking mechanism 
in a non perturbative fashion [23]. Also from some related quark-loop calculations with 
qq mesons in recent times [24,28], it was found that for such 'low energy' processes the 
mass function m(p) is rather well approximated by m(0), so that [22], m q = m(0) + m c . 
Therefore the d—u mass difference is the same at the 'constituent' or at the 'current' levels, 
and this is what has been denoted by Sm in the text (figs. 2). Its smallness compared to m q 
justifies its neglect in all the functions except where it appears explicitly, viz., fig.2. We 
take its value at 5m = 4MeV, as in related calculations [24,28], while the other quantities 
are predetermined from qq [8] and qqq [9] spectroscopy: 

m q = 2Q5MeV; (3 2 (N) = 0.052GeV 2 (4.10) 

so that there are no free parameters in the entire calculation. The results from Sees. 2-4 
are now summarized for (n-p) as : 

N B 2 = 5.5209 x lO^GeV- 10 ; \e { = 1/3] (4.11) 

SM st = +l.nUMeV; 5M 7 = -0.4396Me\/. (4.12) 

Hence 

5M(net) = +1.28MeV; (vs.l.29MeV : Expt) (4.13) 
which is the principal result of this investigation. 

5 Summary and Conclusion 

This calculation fills up an important gap in the two-tier BSE formalism under 3D kernel 
support (CIA) for a simultaneous investigation of spectra and transition amplitudes of 
both qq and qqq varieties within a single unified framework [12,13]. To recapitulate the 
main points, the (first stage) 3D reductions of both the 2-body and 3-body BSE's had 
yielded good agreement with the respective spectra [8,9], with a common set of parameters 
Cq = 0.27 and uj = 158MeV characterizing the structure of the non-perturbative gluon 
propagator, since a third parameter, the 'constituent' mass m q needed for spectroscopy 
[8,9], is essentially the dynamical mass function m(p) in the low momentum limit [22, 23]. 

More substantial tests of the formalism are expected from the (second stage) recon- 
struction of the 4D hadron-quark vertex function which carries the non-perturbative off- 
shell information in a closed form. This exercise was initially confined to the meson-gg 
vertex function whose exact reconstruction [12] had led to several useful checks, from 
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4D loop integrals for hadronic and e.m. transition amplitudes [11,12], to like integrals 
probing the momentum dependence of the quark mass function m{p) which is the 'chiral' 
limit (M n = 0) [17,21,23] of the pion-quark vertex function. Indeed m{p) acts as the form 
factor for loop integrals determining the vacuum to vacuum transitions, and is found to 
predict correctly several condensates, from the basic < qq > [23] to 'induced' condensates 
[39], under one roof. Further tests of the hadron-quark vertex function have come from 
SU(2) breaking effects like p — u mixing [24] and mass splittings in pseudoscalar mesons 
[28], with only one additional parameter representing the d — u mass difference; see also 
[27]. 

The last link in our two-tier formalism has been a reconstruction of the 4D baryon-ggg 
vertex function on the lines of the 2-body case [12], to facilitate the evaluation of corre- 
sponding 3-body loop integrals. This quantity was conjectured some time ago [13], but 
has only recently found a rigorous derivation [29] in terms of Green's function techniques, 
whose main steps are sketched for completeness in Appendix B. The present work repre- 
sents a first application of this quantity, choosing as an example the problem of the n — p 
mass difference. The physics of this problem is two-fold: i) the qqq vertex function which 
is entirely determined by the same dynamics as for the qq case; ii) the strong and e.m. 
SU(2) breaking, (figs. 2 and 3), considered on identical lines to the corresponding problem 
of pseudoscalar SU(2) mass splittings [28]. As such we have refrained here from giving 
fuller references on the d — u mass difference, most of which originated from Weinberg's 
famous paper [40], but several more references on the physics of the problem may be 
found in [28] . On the other hand the entire derivation shows that no free parameters are 
involved, so that the final figure (4.13), although a single number, must not be treated as 
an isolated quantity, but as an integral part of a much bigger package. 

As a point of detail, we should also mention our neglect of the condensate contributions 
inserted in the internal quark lines as in figs. 2, in view of a recent finding [28] that such 
contributions are small within the BSE framework. This may not be too surprising since, 
unlike in the QCD-SR method [26] where such condensate contributions are the principal 
source of non-perturbative effects, this is no longer the case in the present BSE treatment 
which was primarily designed to incorporate non-perturbative effects in the zeroeth order. 
(As a result, the condensate effects in the present non-perturbative scenario may well be 
residual). In this respect the good agreement of our net estimate of 8M, eq.(4.13), with 
experiment (without free parameters) suggests a good support for the BS dynamics when 
viewed together with other related phenomena [8,9,12,23,24,39] within the same dynamical 
framework [12,13,23,29]. 

We also recapitilate the 'QCD' status of this BSE formalism [12] viv-a-vis contem- 
porary methods like QCD sum rules [26] or chiral perturbation theory [27]. As already 
explained in Sec.l, the gluon exchange character of the pairwise qq or qq interactions 
lends them a natural chiral invariance property at the input Lagrangian level with 'cur- 
rent' quarks. In particular, the 'constituent' mass is not an input, but emerges as the 
low momentum limit of the dynamical mass function m(p) that characterizes the quark 
propagators appearing in the 2- and 3- body BSE's, as a result of DBxS [21, 17, 19, 23], 
since the 'current' masses of u, d quarks give only a small additive contribution [22]. The 
empirical aspect of the gluon propagator concerns only its non-perturbative regime which 
often requires separate parametrization even in orthodox formulations [18]. In the present 
formulation, its explicit parametrization with two constants C and ujq) [23] is the price 
for a 'closed form' representation of non-perturbative effects in the derived hadron-quark 
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vertex function, but the returns are rich, e.g., the description of many items that depend 
sensitively on the details of such form factors in the quark-loop integrals [12,23,24,28], with 
their structures firmly rooted in spectroscopy [8,9]. In contrast, the chiral perturbation 
theory [27] has a more explicit QCD content, but with greater emphasis on a perturbative 
treatment, as revealed by expansions in powers of small momenta and "current" masses 
m c [27] for a systemic derivation of the low energy structure of the Green's function in 
QCD [27]. It is a powerful method, highly successful in predicting items like ground state 
masses as well as their splittings, but its lack of a closed form representation prevents 
an equally successful prediction of 'soft' QCD effects in enough details, such as the mo- 
mentum dependence of the mass function, or of hadron-quark vertex functions in general, 
with other observable consequences such as failure to predict the L-excited spectra [10]. 
The method of QCD sum rules [26] also shares similar features such as lack of closed form 
representation for the form factors, and while it does predict items like the n — p mass 
difference, the highly indirect nature of such derivations [41] brings out the parametric 
uncertainties involved in the simulation of soft QCD effects. 

Finally we should like to comment on the principal motivation for this investigation, 
namely to demonstrate the practical feasibility of such realistic quark-loop calculations for 
the relativistic 3-quark problem with a full-fledged (BS) dynamical framework whose basic 
parameters are linked all the way to spectroscopy. The present calculation indeed suggests 
that not only quark-loops involving mesons [12,23,24,28] but even those involving the (less 
trivial) qqq baryon are amenable to a similar degree of dynamical sophistication without 
excessive efforts, so that it makes sense to speak of an effective "4-fermion coupling" 
for both qq and qq pairs within a common parametric framework. This is somewhat 
remiscient of Bethe's "second principle" theory, originally suggested at the two-nucleon 
level of nuclear forces, now reinterpreted at the quark level, with a simple extension to 
include the antiquark in the dynamical description. (This extension would not make sense 
at the NN level, since the NN and NN forces are very different from each other). Indeed 
such a dynamics had been strongly suggested (with concrete examples) in a perspective 
review not too long ago [42], but it seemed to have gone largely by default, as evidenced 
by a strong tendency in the contemporary literature to continue to rely on "ad-hoc form 
factors" [43] to simulate the vertex functions, instead of generating them dynamically. 
Hopefully, some efforts in this direction have been recently in evidence [44], using the 
Nambu Jonalasino model [21] of contact 4-fermion interactions. It is to be hoped that 
Bethe's "second principle" perspective will be upheld by such investigations, at least until 
such time as a fully satisfactory solution to QCD is forthcoming. 

The initial draft of this paper was prepared at the National Institute of Advanced 
Studies. We are grateful to Dr. Raja Ramanna for the warm NIAS hospitality. We also 
acknowledge Ms Chandana's help with some 'difficult' figures. 

Appendix A: Evaluation of the Integral (4.5) 

The master expression (4.8) after being substituted in the full e.m. self energy contribution 
(4.5) is integrated over each dau and doyi ■ The final result is 
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Ftfw, q'i2,P3)exp(-^[q 2 12 + q' 2 2 + 3pt]/f3 2 ) (A.l) 



where 
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(m q + u 12 )(m q + u l2 ')k 2 + (Mw 3 - ^(M 2 - m 2 ))[K, + u 12 ) 2 

+ {m q + cji 2 ') 2 + (m q + wi 2 )(m 9 + uj 12 ) - (m q + u^D^' /2ui 2 ' 
-(m q + u 12 ') 2 D 12 /2ui 2 - (m q + u 12 )(m q + cji 2 ')[£>i 2 /2^i 2 
+D 12 72cu 12 ']A; 2 [(m (/ + 2uo 12 ){m q + 2w 12 ') - m 2 ]/2 
-[(m q + 2cu 12 )(m g + 2u 12 ') - m 2 q ][D 12 /2uo 12 + D 12 ' /2u l2 ']/2 

1 D\ 2 D\ 2 , 4 , 2 



8 6^12^12 



7 + A; - fc 2 [A 2 M 2 + iVM 2 '] (A.2) 



Using eq.(4.9), the integration over q 12 and <^ 2 can be done independently of each other, 
and thus can be written in a compact notation as follows 

SW = £(+^ 2 t 3 ) I f,d^±^L Fl e l-Pim (A.3) 

123 9 J ZUJ 3 

where 

Fi = J n Jn + [Mu; 3 -^M 2 + ^m 2 q ]{2J 20 J 00 + 2J w J w )-2J 20 I 00 
— 2 J10/10 + J11 ' J\i '/2 — Joi<^oi m g/2 — ho'Jio' 

+m 2 q I 00 J 00 - Joo^oo/2 - ^02^02/4 + Joi^oi (A.4) 
and with (n = 0, 1, 2; m = 0, 1, 2), 

J nm ; J nm = 2- 1 / 2 j g 12 rfg 12 e<-§&//3 2 )[V2g 12 p(m, + u 12 ) n [l; ^ 12 /o; 12 ]; (A.5) 



J nm ';/ nm ' = 2" 1 / 2 ! q 12 dq 12 e(-i^ a )[v^g 12 ] m (m, + 2w 12 ) n [l; ^12/^12]; (A.6) 

Appendix B: Derivation of the qqq Vertex Structure, 
Eq.(2.5) 

B.l: Method of Green's Functions 

We outline here some essential steps leading to a formal derivation of eq.(2.5) which was 
written down in a semi-intuitive fashion in [13]. To that end we shall employ the method 
of Green's functions for 2- and 3- particle scattering near the bound state pole, since the 
inhomogeneous terms are not relevant for our purposes. For simplicity we shall consider 
identical spinless bosons, with pairwise BS kernels under CIA conditions [12], first for the 
2-body case for calibration, and then for the 3-body system. 
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B.2: Two-Quark Green's Function 

Apart from some results already goven in the text, we shall use the notation and phase 
conventions of [12,13] for the various quantities (momenta, propagators, etc). The 4D 
qq Green's function G(pip 2 ;pi'p 2 ') near a bound state satisfies a 4D BSE without the 
inhomogeneous term, viz. [12,13], 

i(2n) i G(p 1 p 2 ;p 1 'p 2 ') = Af 1 ^" 1 J dp 1 "dp 2 "K(p 1 p 2 ;p 1 "p 2 ")G(p 1 "p 2 "-p 1 'p 2 ') (B.2.1) 

where 

A 1 =p 1 2 + m q 2 , (B.2.2) 

and m q is the mass of each quark. Now using the relative 4- momentum q= ij>\ — p 2 )/2 
and total 4-momentum P = pi + p 2 (similarly for the other sets), and removing a 5- 
function for overall 4-momentum conservation, from each of the G- and K- functions, 
eq.(B.2.1) reduces to the simpler form 

i(27i) A G(q.q') = Ar 1 A 2 - 1 J dq"Mda"K(q, q")G(q" ', q') (B.2.3) 

where % = q^ — crP^, with a = (q.P)/P 2 , is effectively 3D in content (being orthog- 
onal to Pfj). Here we have incorporated the ansatz of a 3D support for the kernel K 
(independent of a and a'), and broken up the 4D measure dq" arising from (2.1) into the 
product dq"Mda" of a 3D and a ID measure respectively. We have also suppressed the 
4-momentum P^ label, with (P 2 = — M 2 ), in the notation for G(q.q'). 
Now define the fully 3D Green's function G(q, q') as [29] 

G{q,q') = J J M 2 dado-' G{q,q') (B.2.4) 
and two (hybrid) 3D-4D Green's functions G(q,q'), G(q,q') as 

G(q, q') = J MdaG(q, q'); G(q, q') = J Mda'G(q, q'); (B.2.5) 
Next, use (B.2.5) in (B.2.3) to give 

l (2n) 4 G(q, q') = A 1 ~ 1 A 2 ~ 1 J dq"K(q, q")G(q", q') (B.2.6) 
Now integrate both sides of (B.2.3) w.r.t. Mda and use the result [12] 

J MdaA^A^ 1 = 2mD- 1 (q); D(q) = Au(Cj 2 - M 2 /4); u 2 = m 2 + q 2 (B.2. 7) 
to give a 3D BSE w.r.t. the variable q, while keeping the other variable q' in a 4D form: 

(2n fG(q, q') = D' 1 J dq"K(q, q")G(q", q') (B.2.8) 

Now a comparison of (B.2.3) with (B.2.8) gives the desired connection between the full 
4D G-function and the hybrid G(q, g')-function: 

2mG(q, q') = D(q) A^ 1 A^Giq, q>) (B.2.9) 
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Again, the symmetry of the left hand side of (B.2.9) w.r.t. q and q' allows us to write the 
right hand side with the roles of q and q' interchanged. This gives the dual form 

2niG(q, q') = D(g')A 1 , - 1 A 2 , - 1 G(g, q') (B.2.10) 

which on integrating both sides w.r.t. Mda gives 

2mG(q, q') = D^A^A^Giq, q'). (B.2.11) 

Substitution of (B.2.11) in (B.2.9) then gives the symmetrical form 

(2m) 2 G(q, q') = D^Ar'A^Giq, q')D(q') A/" 1 A 2 '^ (B.2.12) 

Finally, integrating both sides of (B.2.8) w.r.t. Mda', we obtain a fully reduced 3D BSE 
for the 3D Green's function: 

(2vr) 3 G(g, q') = D~\q J dq" 'K(q,q")G(q" \q') (B.2.13) 

Eq. (B.2.12) which is valid near the bound state pole (since the inhomogeneous term has 
been dropped for simplicity) expresses the desired connection between the 3D and 4D 
forms of the Green's functions; and eq(B.2.13) is the determining equation for the 3D 
form. A spectral analysis can now be made for either of the 3D or 4D Green's functions 
in the standard manner, viz., 

G(q, q') = £ $ n (g; P)V n {J- P)/(P 2 + M 2 ) (B.2.14) 

n 

where $ is the 4D BS wave function. A similar expansion holds for the 3D G-function 
G in terms of 4> n (q)- Substituting these expansions in (B.2.12), one immediately sees the 
connection between the 3D and 4D wave functions in the form: 

2iri$(q, P) = Ar 1 A 2 - 1 D(g)0(g) (B.2.15) 

whence the BS vertex function becomes T = D x (f>/(2m) as found in [12]. We shall make 
free use of these results, taken as qq subsystems, for our study of the qqq G-functions in 
Sections 3 and 4. 

B.3: 3D Reduction of the BSE for 3-Quark G-function 

As in the two-body case, and in an obvious notation for various 4-momenta (without 
the Greek suffixes), we consider the most general Green's function G(pip2Ps', P1P2P3') for 
3-quark scattering near the bound state pole (for simplicity) which allows us to drop the 
various inhomogeneous terms from the beginning. Again we take out an overall delta 
function 5{p\ + p 2 + p% — P) from the G-function and work with two internal 4-momenta 
for each of the initial and final states defined as follows [13]: 

= pi - pi ; 3?73 = -2p 3 +P1+P2 (B.3.1) 

P = Pl+P2+P3=Pl' +P2 +P3 (B.3.2) 
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and two other sets £1, rji and £2, V2 defined by cyclic permutations from (B.3.1). Further, as 
we shall consider pairwise kernels with 3D support, we define the effectively 3D momenta 
pi, as well as the three (cyclic) sets of internal momenta £i,f)i, (i = 1,2,3) by [13]: 



Pi=Pi- ; ii = i%- SiP ; fji - UP 



(B.3.3) 




(B.3.4) 
(B.3.5) 



Pi + P2 + P3 = ; v x + v 2 + is 3 = 1 



(B.3.6) 



Strictly speaking, in the spirit of covariant instantaneity, we should have taken the relative 
3D momenta £,77 to be in the instantaneous frames of the concerned pairs, i.e., w.r.t. the 
rest frames of Py = pi + Pj] however the difference between the rest frames of P and 
is small and calculable [13], while the use of a common 3-body rest frame (P = 0) lends 
considerable simplicity and elegance to the formalism. 

We may now use the foregoing considerations to write down the BSE for the 6-point 
Green's function in terms of relative momenta, on closely parallel lines to the 2-body case. 
To that end note that the 2-body relative momenta are q^ = (pi — pj)/2 = >/3£fc/2, where 
(ijk) are cyclic permutations of (123). Then for the reduced qqq Green's function, when 
the last interaction was in the (ij) pair, we may use the notation G(£,kVk]£,kVk'), together 
with 'hat' notations on these 4-momenta when the corresponding time-like components 
are integrated out. Further, since the pair £kiVk is permutation invariant as a whole, 
we may choose to drop the index notation from the complete G-function to emphasize 
this symmetry as and when needed. The G-function for the qqq system satisfies, in the 
neighbourhood of the bound state pole, the following (homogeneous) 4D BSE for pairwise 
qq kernels with 3D support: 



i(2n) 4 G(Z m ?rf) = £ AfV / dq'( 2 Mda 12 "K(q 12 , £' 2 )G(& V; 6 V) (B.3.7) 



where we have employed a mixed notation (q 12 versus £ 3 ) to stress the two-body nature 
of the interaction with one spectator at a time, in a normalization directly comparable 
with eq.(B.2.3) for the corresponding two-body problem. Note also the connections 



The next task is to reduce the 4D BSE (B.3.7) to a fully 3D form through a sequence 
of integrations w.r.t. the time-like momenta Sj,tj applied to the different terms on the 
right hand side, provided both variables are simultaneously permuted. We now define the 
following fully 3D as well as mixed (hybrid) 3D-4D G-functions according as one or more 
of the time-like variables are integrated out: 



123 



012 = v^/^; gi2 = V3&/2] fj 3 = -p 3 , etc 



(B.3.8) 




(B.3.9) 



which is S^-symmetric. 




(B.3.10) 
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G^(iv,N) = J Jds 3 ds 3 'G(Zri;Z'ri'y, (B.3.11) 

The last two equations are however not symmetric w.r.t. the permutation group S3, since 
both the variables £, rj are not simultaneously transformed; this fact has been indicated in 
eqs.(B.3.10-ll) by the suffix "3" on the corresponding (hybrid) G-functions, to emphasize 
that the 'asymmetry' is w.r.t. the index "3". We shall term such quantities "S'3-indexed" , 
to distinguish them from S^-symmetric quantities as in eq.(B.3.9). The full 3D BSE for the 
G-function is obtained by integrating out both sides of (B.3.7) w.r.t. the st-pair variables 
dsids/dtidt/ (giving rise to an ^-symmetric quantity), and using (B.3.9) together with 
(B.3.8) as follows: 

(2vr) 3 G(^') = Y.D-\qi2) [ d& 2 K{q 12 ^l 2 )G{^^v) (B.3.12) 

123 

This integral equation for G which is the 3-body counterpart of (B.2.13) for a qq system 
in the neighbourhood of the bound state pole, is the desired 3D BSE for the qqq system in 
a fully connected form, i.e., free from delta functions. Now using a spectral decomposition 
for G 

G(ifj;i'v') = J2MimPWn(i'v;P)/(P 2 + M 2 ) (B.3.13) 

n 

on both sides of (B.3.12) and equating the residues near a given pole P 2 = —M 2 , gives 
the desired equation for the 3D wave function for the bound state in the connected 
form: 

(2tt) 3 0(^;P) = Y J D-\q 12 ) f d^ 2 K{q 12 ^l 2 )<t>{i"f ) "-P) (B.3.14) 

123 

Now the SVsymmetry of in the (£j, fji) pair is a very useful result for both the solution 
of (B.3.14) and for the reconstruction of the 4D BS wave function in terms of the 3D wave 
function (B.3.14), as is done in the subsection below. 

B.4: Reconstruction of the 4D BS Wave Function 

We now attempt to re-express the 4D G-function given by (B.3.7) in terms of the 3D G- 
function given by (B.3.12), as the qqq counterpart of the qq results (B. 2. 12-13). To that 
end we adapt the result (B.2.12) to the hybrid Green's function of the (12) subsystem given 
by G^, eq.(B.3.10), in which the 3-momenta 7)3, 773 play a parametric role reflecting the 
spectator status of quark #3, while the active roles are played by qi 2 , q\ 2 = V3(^3, £3')/% 
for which the analysis of subsec.B.2 applies directly. This gives 

(2«) 2 G 3 ,(^3;6 %) = P(gi 2 )Ar 1 A 2 - 1 G(^ 3 ; i&)Dtf 12 ) A 1 '- 1 A 2 '- 1 (B.4.1) 

where on the right hand side, the 'hatted' G-function has full S'3-symmetry, although (for 
purposes of book-keeping) we have not shown this fact explicitly by deleting the suffix 
'3' from its arguments. A second relation of this kind may be obtained from (B.3.7) by 
noting that the 3 terms on its right hand side may be expressed in terms of the hybrid 
G 3 £ functions vide their definitions (B.3.11), together with the 2-body interconnection 
between (£3,^3') and (^3,^3) expressed once again via (B.4.1), but without the 'hats' on 
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r/3 and 773'. This gives 

(v^) 2 G(&ife;&V) = (yfari) 2 G(tr,;?rf) 

= ^Ar 1 A 2 - 1 (^v / 3) [dq'; 2 Mda 12 ''K(q 12 ,q'( 2 )G(^V^ 3 'rj^ 

123 

= 2 ^(912) Ar 1 A 2 - 1 G 3e (| 3 773; Ai'-'Aa'" 1 (B.4.2) 

123 

where the second form exploits the symmetry between £,77 and 

At this stage, unlike the 2-body case, the reconstruction of the 4D Green's function 
is not yet complete for the 3-body case, as eq. (B.4.2) clearly shows. This is due to the 
truncation of Hilbert space implied in the ansatz of 3D support to the pairwise BSE kernel 
K which, while facilitating a 4D to 3D BSE reduction without extra charge, does not have 
the complete information to permit the reverse transition (3D to 4D) without additional 
assumptions; see [29] for details. The physical reasons for the 3D ansatz for the BSE 
kernel have been discussed in detail elsewhere [23,29], vis-a-vis contemporary approaches. 
Here we look upon this "inverse" problem as a purely mathematical one. 

We must now look for a suitable ansatz for the quantity G 3 £ on the right hand side 
of (B.4.2) in terms of known quantities, so that the reconstructed 4D G-function satisfies 
the 3D equation (B.3.12) exactly, as a "check-point" for the entire exercise. We therefore 
seek a structure of the form 

Gac&ife; i'sVs) = G(Lm; &&) x F(p 3 ,p 3 ') (B.4.3) 

where the unknown function F must involve only the momentum of the spectator quark 
#3. A part of the r) 3 ,r) 3 dependence has been absorbed in the G function on the right, 
so as to satisfy the requirements of S^-symmetry for this 3D quantity [29]. 

As to the remaining factor F, it is necessary to choose its form in a careful manner so 
as to conform to the conservation of 4-momentum for the free propagation of the spectator 
between two neighbouring vertices, consistently with the symmetry between p 3 and p 3 . 
A possible choice consistent with these conditions is the form (see [29] for details): 

F(P3: Ps) = C 3 A 3 - l 5(u 3 - v 3 ') (B.4.4) 

Here A 3 _1 represents the "free" propagation of quark #3 between successive vertices, 
while G 3 represents some residual effects which may at most depend on the 3-momentum 
p 3 , but must satisfy the main constraint that the 3D BSE, (B.3.12), be explicitly satisfied. 

To check the self-consistency of the ansatz (B.4.4), integrate both sides of (B.4.2) w.r.t. 
ds 3 ds 3 'dt 3 dt 3 to recover the 3D ^-invariant G-function on the left hand side. Next, in 
the first form on the right hand side, integrate w.r.t. ds 3 ds 3 on the G-function which 
alone involves these variables. This yields the quantity G 3 ^. At this stage, employ the 
ansatz (B.4.4) to integrate over dt 3 dt 3 . Consistency with the 3D BSE, eq. (B.3.12), now 
demands 

C 3 j J dv 3 dv 3 'A 3 - l 5{v 3 - u 3 ') = 1; (sincedt = dv) (B.4.5) 

The ID integration w.r.t. dv 3 may be evaluated as a contour integral over the propagator 
A -1 , which gives the pole at v 3 = u> 3 /M, (see below for its definition). Evaluating the 
residue then gives 

G 3 = »/(Mw 3 ); Col = m q 2 + p 2 3 (B.4.6) 
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which will reproduce the 3D BSE, eq.(B.3.12), exactly] Substitution of (B.4.4) in the 
second form of (B.4.2) finally gives the desired 3-body generalization of (B.2.12) in the 
form 

3G(^;eV) = E^(^2)AiFA 2F D(gi 2 )A 1F 'A 2F 'G(^ 3 ;^3)[A3F/(M7r^3)] (B.4.7) 

123 

where for each index, = —iA^ 1 is the Feynman propagator. 

To find the effect of the ansatz (B.4.4) on the 4D BS wave function P), we do 

a spectral reduction like (B.3.13) for the 4D Green's function G on the left hand side of 
(B.4.2). Equating the residues on both sides gives the desired 4D-3D connection between 
$ and 0: 

From (B.4.8) and eq.(2.1) of the text, we infer the structure of the baryon-ggg vertex 
function V3 as given in eq.(2.5) of the text. For a detailed discussion of the significance 
of this result, vis-a-vis contemporary approaches, see [29]. 

Figure Captions 

Fig.l: Diagrams for BS normalization of Baryon-ggg vertex function. 1(a) shows quark 
#1 emitting a zero momentum photon (k — 0); its last qq interaction was with #2, while 
#3 is the spectator. 1(b) is the same diagram with the roles of #1 and #2 interchanged. 
1(c) denotes schematically two more such pairs of diagrams obtained with cyclical per- 
mutations of the indices (123) in pairs. The 4- momenta on the quark lines are shown as 
used in the text. 

Fig. 2: Diagrams for the two-point interactions of the quark lines with the mass shift 
operator —8mT^/2 in place of the photon in fig.l, but otherwise with identical topolog- 
ical correspondence of figs.2(a,b,c) to figs.l(a,b,c). 

Fig. 3: Diagrams for the e.m. self-energy of the uud (proton) configuration. 3(111) is 
shown in detail with full momentum markings as employed in the text, and corresponds 
to quark #3 as the spectator, while the quark lines #1 and #2 are joined by a transverse 
photon line. Similarly 3(1) and 3(11) correspond to #1 and #2 respectively as spectators 
in turn. Note that, unlike in fig.l and fig. 2, the interchange of #1 and #2 in fig.3(III) 
does not give a new configuration. 

Note : Due to lack of adequate software for a proper 'DVT rendering of the three 
figures above, it has not been found possible to include the same in this file. The inconve- 
nience on this account is regretted. However, any interested reader will find it quite easy 
to reconstruct the three "quark loop" figures (in each of which there are three internal 
quark lines between two baryon vertex blobs) on the basis of the information supplied in 
the captions above. As a further guidance, the interested reader may refer to ref. [28] 
which contains similar figures with two quark lines each between the corresponding vertex 
blobs of the pseudoscalar mesons concerned. 
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